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The recently proposed Partition Theory (PT) (J. Phys. Chem. A 2007, 111, 2229), is illustrated 
on a simple one-dimensional model of a heteronuclear diatomic molecule. It is shown that a sharp 
definition for the charge of molecular fragments emerges from PT, and that the ensuing population 
analysis can be used to study how charge redistributes during dissociation and the implications of 
that redistribution for the dipole moment. Interpreting small differences between the isolated parts' 
ionization potentials as due to environmental inhomogeneities, we gain insight into how electron 
localization takes place in H J as the molecule dissociates. Furthermore, by studying the preservation 
of the shapes of the parts as different parameters of the model are varied, we address the issue of 
transferability of the parts. We find good transferability within the chemically meaningful parameter 
regime, raising hopes that PT will prove useful in chemical applications. 



1. INTRODUCTION 

Consider a molecule of composition AB with parts A 
and B having different ionization potentials when iso- 
lated. A long-standing problem is how to associate 
charges with each part as the parts are separated. At 
intermediate separations, one expects that, as bonding 
electrons would spend unequal time in the vicinity of each 
part, one would have to assign non-integer average num- 
bers of electrons to each, numbers which become integers 
at infinite separation. Density-functional theory (DFT) 
is defined only for integer electron numbers as originally 
developed [J @j- If the dependences of the energy func- 
tionals of integer DFT on electron density were contin- 
ued to densities containing noninteger electron numbers 
and applied to the separation of AB into A + B, at in- 
finite separation A and B would have unphysical nonin- 
teger electron numbers, as pointed out by Perdew et. al. 
(PPLB) 0,11]. Instead, PPLB argued that an ensemble 
generalization of ground-state DFT should be used for 
systems with noninteger electron number. 

With how to treat noninteger electron number resolved 
by PPLB, the issue of how rigorously and systematically 
to decompose a system into its parts remains [5, 6]. Two 
of the present authors have proposed an exact scheme, 
partition theory (PT) @, H, 0], based on the PPLB en- 
semble DFT. In Q, their PT was brought to full formal 
development and used for a reconstruction of chemical 
reactivity theory which eliminated the inconsistencies of 
earlier formulations and enriched them. Applying PT to 



the case introduced above, AB — > A + B, the parts would 
obviously be A and B. 

To illustrate the conceptual structure and physical con- 
tent of PT, a very simple system was studied in [Io| . 
a caricature of the hydrogen molecule consisting of two 
electrons moving independently in one dimension under 
the influence of two attractive delta-function potentials 
of equal strength, lDi?2- In the present paper, a cor- 
responding model of a heteronuclear diatomic molecule 
AB is studied via PT (1DAB). The model once again 
consists of two electrons moving independently in one 
dimension under the influence of two attractive delta- 
function potentials of unequal strengths — Za and — Zb 
with Zb < Za- Since a one-electron _B-atom would tend 
to donate its electron to the more electronegative ^4-atom 
when brought together, A can be thought of as a Lewis 
acid and B as a Lewis base. 

In the limit Za = Zb = Z, the model becomes that 
treated in pjj • Moreover, reducing the number of elec- 
trons from 2 to 1 requires little modification of that an- 
alytic theory, and its numerical results can be used to 
examine the dependence of charge transfer on Za — Zb 
and inter "nuclear" separation. The resulting theory can 
also be used to explore how symmetry breaking local- 
izes the single electron of when it is separated into A 
and A + , a subtler problem than the localization of both 
electrons on A when AB — ► A~ + B + . 

In Section 2, the formalism developed for the IDH2 
problem in ref. [T(| is extended to the present 1DAB prob- 
lem. Numerical results are given in Section 3 for the de- 
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pendence of the electron densities of the parts, for the 
charge transfer from A to B, and for the dipolc moment 
as functions of Za, Zb, and the inter "nuclear" separation 
R. Results are also given there for the partition potential 
of PT, and how it induces electronegativity equalization 
between the parts is discussed. The transferability of the 
properties of the atoms is discussed as well. The united- 
atom limit, of more academic interest than chemical rele- 
vancy, is discussed separately at the end of Section 3. Sec- 
tion 4 is devoted to the one-electron molecule ID AB + . 
The case Za I Zb is used to show how trivial symme- 
try breaking, {Za — Zb) « Za, is sufficient to localize 
the electron on A, illustrating how real separates 
into H + H + because of small environmental perturba- 
tions. We conclude in Section 5 with a brief discussion of 
the significance of these very simple illustrations of the 
power and utility of PT for population analysis and for 
the transferability of fragments with their properties be- 
tween different molecular contexts. Detailed derivations 
of all analytic results are presented in an Appendix. 



2. 1DAB; INDEPENDENT ELECTRONS 
MOVING IN UNEQUAL 5-FUNCTION 
POTENTIALS IN ONE DIMENSION 

A. The Molecule 

In we considered an analogue of the Hi molecule 
in which two electrons move independently in one dimen- 
sion under the influence of two ^-function potentials of 
equal strength —Z. In the present section, we consider 
the heteronuclear analogue 1DAB in which the nuclear 
(^-function of the acid A is of strength —Za and that 
of the base of strength —Zb, with Za > Zb- These 
"nuclear charges" are allowed to vary continuously. The 
ground-state wave functions ipa( x ) an d energies E^ of 
the isolated "atoms" a — A,B are (atomic units are used 
throughout): 

iPlix) = ^/Z~exp[-Z a \x\] , (2.1) 

E° a = -Zl/2 . (2.2) 

The ground-state energy Em{Nm = 1) of one elec- 
tron moving independently in the two <5-functions, that 
of strength —Za at x = —R/2 and that of strength —Zb 
at x = R/2, is —k 2 /2, where 

{k- Z A ){n - Z B ) = e- 2KR Z A Z B ■ (2.3) 

The solutions of Eq. (|2.3p are plotted as a function of 
internuclear separation in Subsection A of the Appendix. 
From here on we are concerned only with the lowest- 
energy solution, denoted simply as k, which corresponds 
to a bonding state that is doubly occupied when Nm = 2. 




FIG. 1: Molecular density for three different internuclear sep- 
arations: R = 0.4 (left), R = 2.0 (center), and R = 4.0 
(right). In this plot, Z A = 1.02, and Z B = 0.98. 



The corresponding ground-state wave function, ipM (x) 

is 

$ M {x) = Ce< R / 2+x ^ , x<-R/2 , 

= De KX + Fe~ KX , - R/2 < x < R/2 (2.4) 
= Ge K{R / 2 - x ^ , R/2<x , 

where C , D, F, and G are constants whose explicit ex- 
pressions in terms of Z a, Zb, R, and n are given in 
subsection A of the Appendix fEcis lA.ll IA.2[) . In the 
combined-atom limit R J, 0, G = C = k 1 / 2 , and D 
and F are irrelevant. In the limit R f oo, C — k 1 / 2 , 
D = e-^/V/ 2 , and G = F = so that ip M {x) is local- 
ized on A only. 

The two-electron molecular density is 

n M (x) = 2\4> M {x)\ 2 ; (2.5) 

the total energy of the molecule M is 

E M (N M = 2) = 2E M (N M = 1) = -k 2 ; (2.6) 

and the chemical potential of the molecule is 

Hm=E m {2)-E m (1) = -k 2 /2 , (2.7) 

all just as for the lD_ff2 case of ref. [10]. FigfT] shows the 
dependence of um{x) on internuclear distance R for Za — 
1.02, Zb = 0.98. For these only slightly different values of 
Za and Zb the transition from primarily ionic character 
(right panel) to mixed ionic-covalent character (center 
panel) occurs at relatively large bond length R ~ 3.2. 
The left panel of FigQ] shows a density belonging to the 
interesting but less chemically meaningful united-atom 
regime, a regime that we discuss separately in Subsection 
2D because it allows us to draw conclusions regarding the 
limits of utility of Partition Theory. 

B. The parts 

Our task is to partition um (x) into contributions from 
the two parts of M, fragments A and B, 

n M (x) = n A (x) +n B (x) , (2.8) 

with ua,b localized primarily around —R/2, R/2, respec- 
tively. Because ipM and therefore tim is larger near A 
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than near B (recall that Za > Zb) 
bers of the fragments, 



the electron num- 



dxnA,B(x) 



are unequal with Na > Nb and 

N A + N B = N M = 2 



(2.9) 



(2.10) 



Na,b arc nonintcgcrs in general so that Eq. (|2.1[) implies 
that 



N A 



N F 



V , 0<V<1 



(2.11) 



and the PPLB ensemble [3] must be used for the frag- 
ments. In this use of the PPLB ensemble lies the main 
difference with Parr's atoms-in- molecules approach based 
on a minimum promotion energy criterion |ll| . For A, a 
singly-occupied state of either spin occurs with probabil- 
ity v/2, and a doubly-occupied state occurs with proba- 
bility 1 — v in the PPLB ensemble. For B, an unoccupied 
state occurs with probability 1 — v and a singly-occupied 
state of either spin occurs with probability u/2. The 
densities of the fragments are, accordingly, 

n A {x) = {2-v)^ A {x) , n B {x) = v^ 2 B {x) , (2.12) 

where , >p a (x) (a = A, B) is a real one-electron wavefunc- 
tion localized around —R/2, a = A, or +R/2, a = B. 
Analytic expressions for the rp a (x) are derived in subsec- 
tions B through D of the Appendix. They can be written 
most simply as: 

ip A (x) = (2 - v)- 1/2 ij M (x) (cosp(x) + sin/3(x)) , 

(2.13) 

ip B {x) = v- 1/2 ij) M (x) (cosP(x) - sm/3(x)) , (2.14) 

where (3{x) is an auxiliary polar-angle function that de- 
termines v according to 



1 



dx^ 2 M (x) sin 2[3{x) 



(2.15) 



as shown in detail in Subsection B of the Appendix. The 
corresponding measure of electron transfer q is 



1 



dx^ 2 M (x) sin 2/3(x) 



(2.16) 



In Partition Theory, the wave functions and densities 
of the parts are found by minimizing the sum of the en- 
ergies of the individual parts subject to the constraints 
that the electron densities and numbers of the assembly 
of parts be identical to those of the molecule M. The 
expressions for iPa(x) and 4>b(x) of Eqs. (|2.13p - (|2.14p to- 
gether with Eqs. (|2.1ip - (|2.12p , guarantee that the con- 
straints are met for 1DAB. Only the determination of 
(3{x) remains, and the procedure for determining it by 
minimizing the sum of the energies of the fragments is 
described in Subsections C and D of the Appendix. 




FIG. 2: Molecular density um(x) (solid), and fragment den- 
sities ua(x) (dotted) and ns(a;) (dashed) for Za = 1.05, 
Zb = 0.95 and R — 2. For these values, v = 0.6, indicating 
substantial charge transfer even with only a 10% difference 
between Za and Zb and a small R. 



Numerical evaluation of the quantities of interest, v, 
P{x), and ip a (x), is straightforward for representative Z a 
and various R. For example, Fig. [5]shows the fragments' 
atomic densities for Za — 1.05, Zb = 0.95, and R = 2. 
Interestingly, in spite of having high electron density in 
the bonding region around the molecule's center of mass, 
the fragment densities resemble true atomic densities. 
This is one of the most significant results of our work. 
In order to investigate the extent to which this is true, 
and quantify it, results for different parameter regimes 
are reported and discussed in the next Section. First, 
however, we discuss in the following subsection how the 
partition potential vp(x) of PT can be extracted from 
these results. 



C. The partition potential 

Within the framework of PT [9], the wavefunctions 
ip a {x) for the parts a — A, B satisfy the Schrodingcr 
equation 

{H a + v P {x)) ip a (x) = HMifra(x) , a = A,B , (2.17) 

where the Hamiltonians H a for the parts are defined by 
Eq. (|A.lip . The partition potential vp(x) is the same 
for both parts, as is the eigenvalue [iu corresponding to 
the molecular chemical potential, Eq. (|2.7p . We construct 
vp(x) in subsection E of the Appendix, obtaining: 



\Z 2 A ^ M {-R/2)cos 2 213, 



6{R/2- |x|)(2.18) 



- ^2v a (1 - s a tan/3 Q (l + cos2/3 Q )) , 



where j3 a (a — A, B) are constants that determine the 
upper and lower bounds of f3(x) (derived in subsection D 
of the Appendix), s a = ±1 (plus sign for part A, minus 
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sign for part B), and 9(y)(= for y < 0, 1 for y > 0) is 
the Heaviside step function. Numerical results for vp(x) 
are presented and discussed in the next Section. 



3. ILLUSTRATIVE RESULTS 

For purposes of discussion, we define three critical pa- 
rameters i?ip, i?don, and i? cc as the values of R at which 
significant changes in the ionization potential, the den- 
sity, and the occupation of a fragment, respectively, take 
place. Our definitions for i?ip and i?d cn are immediately 
applicable to any diatomic molecule, and our definition 
for R occ could be easily generalized to be applicable to 
any diatomic molecule. We start by defining R occ in 
Subsection A in order to discuss population analysis and 
charge transfer, and continue in Subsection B with i?den 
and Rip to address the issue of transferability 



A. Population Analysis; Charge Transfer 

First, note that q = 1 — v of an electron is transferred 
from B to A. Figure [3] shows how v changes when R is 
varied for 4 different values of AZ (and Z = 1) for R > 1, 
and Fig [M] does for smaller R as well. Define R occ as the 
larger of the two values of R for which Na = 2Nb (and 
therefore v = N/3), a reasonable criterion for the ionic to 
mixed ionic-covalent crossover. The solid curve in Fig [5] 
displays the behavior of i? occ vs. AZ. We observe that 
i? occ is a sensitive function of AZ , especially as AZ — > 0. 
(We come back to this point in the next paragraph and in 
Section 4 when discussing electron localization in as 
the molecule dissociates) . This is a simple illustration of 
the utility of PT to quantify the degree of charge trans- 
fer taking place as a heteronuclear diatomic molecule is 
stretched out. That charge transfer occurs is readily seen 
from the sequence of molecular densities plotted in FigfTJ 
and PT provides an unambiguous way to characterize it. 

This independence of v is interpretable as a sequence 
of crossovers in the electronic structure of the molecule 
M. At large R, M is ionic with v J. and q f 1. As R 
decreases, there is a crossover around i? OC c from that ionic 
to a mixed ionic-covalent state. As R decreases further 
the covalent character of that mixed state increases. At 
still smaller R, there is another crossover to a combined- 
atom state around the maximum in v. Finally, at very 
small R, there is a rapid crossover back to the ionic state, 
an interesting fact that we discuss further in Subsection 
E. 

It is also of interest to examine how well the prod- 
uct of transferred charge with the internuclear separa- 
tion agrees with the actual electronic component of the 
dipole moment of the molecule. Figure [4] compares this 
point-charge dipole (PC) 

d PC = -qR , (3.1) 
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FIG. 3: Population of fragment B (=u) for various values of 
AZ, for Z = l. 



with the actual electronic part of the dipole moment 
d = J dx xum{x) for various values of AZ. If the 
fragment densities were inversion-symmetric about their 
"nuclei" , there would be exact agreement between dpc 
and d. But there is in addition a dipole moment on 
each fragment from the polarization of its electron cloud, 
d a = J dx xn a (x), d — dpc + d^ + dg. Figure \E\ shows 
d(R) along with the fragment dipoles d a (R), a — A,B, 
for fixed AZ. Clearly, when AZ = 0.04 the fragment 
dipoles are small compared to the molecular dipole only 
for separations larger than R ~ 3. The estimate (|3.1[) 
shows excellent agreement whenever the dipole is signif- 
icant, and only fails (error > 20%) when the dipole is 
small (Fig|5]) . Note also from the inset of Figure H that 
the ^-dependence of the percent difference changes little 
when AZ changes from 0.01 to 0.12. Figure [5] shows that 
dA is of opposite sign to ds and dpc- This sign reversal 
arises from the distortion of ua by covalency, which in- 
creases it in the bonding region. However, in 1DAB dA 
is not large enough in magnitude to overcome ds + rfpc 
even at small R. If it were, d would be of sign opposite 
to that expected from the electronegativity difference of 
A and B. In real diatomic molecules containing small 
bases, such sign reversals are observed [lH, and parti- 
tion theory promises a simple explanation. 



B. Transferability 

The mixed covalent to ionic crossover occurs at quite 
small R for very small AZ when Z = 1. For the crossover 
to occur at an internuclear separation of about lA, AZ 
need only be about 0.09. The shape of each atom, how- 
ever, is not as sensitive. Figure [S] shows iPa(x) (solid) and 
iPb(x) (dotted) for two different separations and three 
different values of AZ covering the same range as that 
of Figj3] It is apparent that orbitals corresponding to 
different values of AZ start differing significantly only 
close to the region where ipA{%) and tpsix) overlap (the 
"bonding" region), and only for large internuclear sep- 



FIG. 4: Electronic dipoles as a function of inter "nuclear" sep- 
aration, both from the point-charge value of Eq. (|3.1|l (solid 
lines), and exactly (dashed), for 3 different values of AZ. The 
inset shows the absolute percent difference. (Z = 1). 




FIG. 5: Electronic dipole moment for the 1DAB molecule 
(solid), and fragment dipole moments, for fixed AZ = 0.04 
and 2 = 1. 



arations R > 3. Since at large intcrnuclear separations 
the orbitals ip a (x) become equal to the isolated solutions 
ipa( x ) °f Eg . i|2. ip . it is interesting to examine how and 
when the ip a (x) start departing from the ip^x). Fig- 
ure [7] shows that for AZ — 0.08, i^a{x) is almost iden- 
tical to ?p A (x) when R = 4.8, differing only slightly in 
its right-hand tail, and that the shape is still preserved 
for chemically-relevant values like R ~ 1.6. In order to 
appreciate the differences more clearly, Fig |8] displays the 
differences of the squares D(x) = ^a(x) 2 — (ip A (x)) 2 . We 
observe that the two orbitals depart appreciably when 
the spatial integral of the absolute value of this quantity 
reaches a value of J dx\D(x)\ ~ 0.2, so we define i?den 
as the corresponding separation. As shown in the inset 
of Fig[71 i?den remains almost constant at i?den ~ 0.66, 
small compared to i? cc in the corresponding range of 
AZ. A similar behavior is observed for i?ip, which we 
define as the value of R at which the ionization potential 
of the molecule begins to differ significantly [20%] from 
the ionization potential of the most electronegative (iso- 



FIG. 6: Fragment orbitals tpA(x) (solid) and iPb(x) (dotted) 
for different values of R (1 and 3). The fragment- A orbitals 
are plotted for three different values of AZ (0.04, 0.08, and 
0.12). Orbitals corresponding to different values of AZ start 
to differ in the overlapping region only for large internuclear 
separations R > 3. Differences due to changes in AZ are 
indistinguishable to the eye for R < 3. 



lated) atom. Figure [5] compares i? occ , Rip and i?don- The 
preservation of the shape of the orbitals in a range of AZ 
where there is signficant charge transfer (i?dcn << -Rocc) 
is a strong indication of the transferability of the frag- 
ments emerging from Partition Theory [13j . The feature 
of transferability can also be seen for fixed AZ by com- 
paring fragments corresponding to different values of R. 
Figure [TU] shows such a comparison for AZ = 0.06 by 
overlapping the iPa(x)'s corresponding to different i?'s 
(this requires shifting the origin of ^4-atoms to the left 
and -B-atoms to the right). We note that even though 
the size of the fragments (controlled by v) changes sub- 
stantially as R varies from 1 to 4, their shape is quite 
insensitive to relatively large changes in R with concomi- 
tantly large changes in um(x) and fiM- For the value of 
AZ chosen in FigfTUl v is only about 0.2 when R = 4, 
and yet the shape of the corresponding fragments only 
differ slightly from those obtained when R — 1 and v is 
close to 0.9. 



C. Partition Potentials 

The corresponding partition potentials are shown in 
Fig. QT] Both parts A and B must have equal electroneg- 
ativities, sharing the same HOMO eigenvalue, Eq. (|2.17D . 
which must be equal to the overall chemical potential 
fiM- The partition potential ensures that this happens 
by acquiring a specific form, with an asymmetric neg- 
ative well in between the fragments and two negative 
delta-functions at ±i?/2. In the limit of infinite separa- 
tion, when k — > Zai the external potential of "nucleus" 
A does not require any correction in order to reach /im, 
but the external potential of "nucleus" B does. Accord- 
ingly, the delta-function component of vp vanishes on A 
at infinite separation, but not on B (fig ll3p . It is negative 
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FIG. 7: Comparison of fragment orbital tpA(x) (solid line) 
and isolated atomic orbital iI>a{ x ) (dotted) for R — 4.8 and 
R = 1.6. The inset shows the slow variation of i?den with AZ 
(see text). 
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FIG. 8: Difference between the fragment density on A and 
the isolated atomic density, (4>a( x )) 2 ~ ^aM) f° r different 
values of R and fixed values of AZ = 0.08 and 2 = 1. 
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FIG. 9: Comparison of i? OC c (value of R at which the occu- 
pation on A is twice the occupation on B), Rip (at which 
the ionization potential of the molecule is 20% larger than 
the ionization potential of an isolated A atom), and i?dcn (at 
which the fragment densities change significantly as compared 
to the corresponding isolated-atom densities - see text). All 
curves with Z = 1. 




FIG. 10: Fragment orbitals tpA(x) and ipB(x) for R = 4 
(solid), R = 2 (dotted), and R = 1 (dashed), when Z A = 1.03 
and Zb = 0.97. The A-fragment corresponding to R = 2 was 
shifted to the left by la.u. and the B-fragment was shifted to 
the right by la.u. The A-fragment corresponding to R — 1 
was shifted to the left by 1.5a.u., and the B-fragment was 
shifted to the right by 1.5a.u. 
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FIG. 11: Left panels: A and B fragments: ipA(x) = 
\J ua(x)/(2 — v) (solid) and i/jb{x) = \fnB~{x)Jv (dashed) 
for R = 4 (upper), R — 2 (center), and R — 1 (bottom). 
Right panels: corresponding partition potentials (^-functions 
at ±R/2 not shown). For all these plots, Za = 1-03 and 
Z B = 0.97. 



there, and has a magnitude smaller than AZ because the 
smooth negative well persists at large separations, getting 
further from A, and closer to B (Figs fTTTfT2|) . thereby con- 
tributing to the lowering of the eigen-energy of part B 
towards hm- Figure [121 shows a closer view of the de- 
pendence of the smooth part of the partition potential 
with AZ and with R, and Figure [T31 illustrates the same 
dependences for the amplitudes of its delta-function com- 
ponents. 
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FIG. 12: Smooth part of the partition potential for 2 different 
values of R and AZ = 0.04 (solid), AZ = 0.08 (dashed) and 
AZ = 0.12 (dotted). For all these curves, Z = 1. 



FIG. 14: Same as Fig[3]with added small- R range to illustrate 
striking united-atom behavior of v (see text). 




FIG. 13: Delta-function weights of the partition potential, 
for AZ = 0.04 (solid), AZ = 0.08 (dashed), and AZ = 0.12 
(dotted). The 3 upper curves correspond to atom A, and the 
3 bottom ones to atom B. For all these curves, 2=1. 
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FIG. 15: Behavior of R' occ for ID AB + when Z = 1. The inset 
shows -Rocc vs - In AZ. The line corresponds to the best linear 
fit in the range —9 < InAZ < —5 (slope: -0.185; intercept: 
-0.502, agreeing with the tight-binding formula derived in the 
Appendix, Eq.((X43j) 



D. United-atom limit 

In order to focus attention on the chemically relevant 
range of internuclear separations, we started Fig[3] at 
R = 1 . The small- R behavior of v can be seen in FigfT4l 
Perhaps counterintuitively, as R decreases below a given 
(small) value, a rapid crossover takes place from a mixed 
ionic-covalent state to an essentially ionic state. In the 
IDH2 model of ref . [lj3] , this reversion to the ionic state is 
absent, the sequence of crossovers being atomic to cova- 
lent to combined-atom state. There is a singularity in the 
1DAB model at AZ = 0, R = 0. If lim fl , A zio {^f } 4 0, 
the combined-atom state persists to R = 0. If, however, 
lim.K_A.zio { ~£z} I 0; a reversion to the ionic state oc- 
curs. This counterintuitive feature illustrates an impor- 
tant limitation to the utility of partition theory. When 
kR becomes significantly less than unity as R decreases, 
the il>a{x) overlap substantially, and the primary moti- 
vation of partition theory, decomposition of the electron 
density into distinct localized components, is frustrated. 



The shape of the fragment densities also departs sig- 
nificantly in this limit from pure atomic densities. For 
example, for the same value of AZ = 0.08 used in FigJSl 
the maximum value of the difference between ipA(x) 2 and 
-0^4 (x) 2 when R = 0.4 is about 6 times larger than the 
corresponding maximum when R — 1.6. 

4. ELECTRON LOCALIZATION IN 
DISSOCIATING IL+ 

Even though the ground-state wavefunction of is 
symmetric, with 50% of its amplitude on the right atom, 
and 50% on the left, the slightest asymmetry due to en- 
vironmental perturbations forces the electron to localize 
onto one of the two nuclei as the molecule dissociates. 
Since this symmetry breaking can now be studied ex- 
perimentally via intense few-cycle laser pulses with con- 
trolled field evolution [l4|], there is resurgent interest in 
theoretical models to describe electron localization dur- 
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ing molecular dissociation (see for example ref.[15[ for 
a recent study of dissociation and ionization of small 
molecules steered by external noise). Our simple the- 
ory of the preceeding sections can be used as such a 
model, provided we interpret small differences between 
the magnitudes of Z A and Zb as due to the effect of 
an inhomogeneous environment. In fact, since we have 
dealt with 2 non-interacting electrons, only minor modi- 
fications of our results are needed in order to analyze the 
one-electron case, lD^li?" 1 " [1(1 : the chemical potential is 
still identical to that given by Eq. (|2.7p . The number con- 
straint of Eq. (|2.10[) is modified to N' M = 1 (we use primed 
symbols to represent 1-clcctron quantities to distinguish 
them from their 2-electron analogs), and Eq. (j2. 1 1[ ) goes 
to: 



1 



N' B 



(4.1) 



The densities of the atoms are: 



n' A (x) = (l-v')ip A (x) , ri B (x) = v'^\(x) , (4.2) 
and following the same steps leading to Eq. (|2.15[) , we find 



v = — 
2 



1 - / dxilj 2 M {x)sm2P'(x) 



(4.3) 



Since the Euler equation for j3(x), Eq. (|A.16p . as well as 
the boundary conditions, remain unchanged, f3'(x) = 
(3(x) and the values that v' takes as a function of R 
and AZ are simply half of those calculated for 1DAB, 
v' = v/2. As before, we define R' occ as the value of R 
for which N' A = 2N' B , corresponding to v' = 1/3 (from 
Eq.gT|) and N' A + N' B = 1). The behavior of R' occ as 
a function of AZ is then identical to that of R occ (AZ). 
The inset of Fig. (TT5|) shows that for small AZ, R occ is 
proportional to In AZ. We conclude that the more inho- 
mogeneous the environment (the larger AZ), the earlier 
electron localization occurs along the dissociation path- 
way, with a logarithmic dependence in this case. As i? cc 
becomes very large when AZ becomes small, this strik- 
ing relation between i? cc and AZ follows from a simple 
tight-binding argument outlined in the Appendix (Sub- 
section F). 



5. CONCLUDING REMARKS 

We have applied Partition Theory Q to a simple model 
of a heteronuclear diatomic molecule. We found ana- 
lytic expressions for the densities of the parts, the charge 
associated with each of the molecular fragments, and 
the partition potential that guarantees electronegativity 
equalization. Numerical calculations for various parame- 
ter regimes allow us to reach important conclusions: (1) 
i?occ has a strikingly different bahavior than Rip and 
i?dcn as AZ — > (FigJS]). Since R occ measures the value 
of inter "nuclear" separation at which significant charge 
transfer occurs, and i?dcn and i?ip measure the value of 



R at which signifcant change in the shape of the frag- 
ment wavefunctions takes place, we conclude that the 
fragments of Partition Theory, at least within this sim- 
ple model, are to large extent transferable. (2) Environ- 
mental fluctuations (modeled by small finite AZ) local- 
ize the single electron of onto one of the two nuclei. 
As dissociates, the more inhomogeneous the environ- 
ment, the earlier localization occurs along the dissocia- 
tion pathway. The explicit results reported both here 
and in ref. [l(| both illustrate important features of Par- 
tition Theory [9[ and support the proposed use of PT for 
a broad range of applications to real systems including 
the sharp definitions of parts of a larger system, popu- 
lation and charge-transfer analysis, and the examination 
of transferability. 



APPENDIX 



A. The Molecule 

Eq. (|2.3p has two solutions for all R < (Za + 
Z B )/2Z A Z B = R*: one, k+, belonging to a bonding 
state doubly occupied when N = 2, and another one, 
«o = 0, belonging to a state at the bottom of the con- 
tinuum. For all R > R* there is another solution, 
corresponding to an unoccupied antibonding state. As 
R | oo, k+ and k_ give rise to the two energies E° A B 
of Eq. (|2.2[) corresponding to the two states ip A B (x) of 
Eq. (|2.ip of which the lower state i(j a (x) localized to the 
Lewis acid A is doubly occupied, the Lewis base B having 
donated its electron. As R J. 0, k + approaches (Za + Zb) 
and gives rise to a doubly occupied combined-atom state 
with the Z a of Eqs.([2TT]) and ([272]) replaced by Z A + Z B . 
The value of k+ decreases monotonically from Za+Zb at 
R I with finite derivative at zero and vanishing deriva- 
tive at infinity. k_ increases monotonically from zero at 
R J, R* to Zb at R t oo. kq remains zero throughout (see 
FigfTSj) . In the text following Eq. (|2.3p . we referred only 
to k + , dropping the subscript + for notational simplicity. 



The i?-dependences of all three solutions of Eq. 
for k are shown in Fig [11] for Za = 1,Z B = 0.9. Note 
that k+ begins to differ significantly from k + (oo) only 
at separation R < R\p ~ 1.8, less than relevant for real 
molecules, except for H2. For example, for Za = 1 and 
Z B = 0.9, one gets from Eqs. (|2.3p and (|2.7p an ionization 
energy equal to that of lithium hydride (LiH), for which 
the equilibrium bond distance is Rq — 3.05 (>> Rip)- 

Finally, the constants of Eq. (|2.4p are given in terms of 
Za, Z b , R, and k, by: 



D 
F 
G 



e KR/2 (1 _ Za/k) Q 

e KR ' 2 {Z a /k) C 
e KR [(« - Z a )/Z b ] C 



(A.l) 
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FIG. 16: k versus _R for with Z A = 1 and -Zg = 0.9. The 
upper curve shows the dependence of k for the doubly occu- 
pied orbital of the ground state which goes over from ionic 
on A at infinite R through mixed ionic-covalent to combined- 
atom at small R. The lower curve shows that of k for the 
empty orbital of the excited state which changes from atomic 
on B at infinite R to antibonding at intermediate R and dis- 
appears for R < R*; ko = is also a solution of Eq. (|2.3[) . The 
vertical dotted line indicates the separation Rip at which k, 
begins to differ significantly (by 10%) from k+(oo). 



XA,b{x) = — ^ipMix) (cos0(x) ± sin/?(a;)) . (A. 9) 

Because Xa,b{x) are non-negative, |/3| cannot exceed 
tt/4. 

Inserting Eq. (|A.9[) into Eq. (|A.5[) and subtracting the 
squares of the two resulting equations leads to an expres- 
sion for v, Eq. (|2.15[) : 

v = 1 - J dayiplf (x) sin 20(x) , (A. 10) 

after integrating over x. Determination of the polar 
angle (3(x) is thus sufficient for the determination of v 
(charge transfer) and the electron population of the frag- 
ments (population analysis). For the symmetric IDH2 
case, P(x) is odd and tpM{x) even so that the integral 
in Eq. (|2.15p vanishes, yielding v = 1 and equally pop- 
ulated fragments. In the present case, since iPm{x) is 
normalized to unity and | sin2/3(x)| < 1 with the domain 
of positive (3(x) weighted more heavily than that of neg- 
ative (3{x), the integral in Eq. (12.15p lies in (0, 1), as must 
v in accordance with Eq. (|2.11[) . 



- e -^ + 2 -^[l + 2R{ K -Z A )]\ 

K K J 

B. The polar angle [3(x) 

Substituting Eqs. l|2TT2|) and [[23]) into Eq. (|2~g|) leads to 
{2-v)i^(x)+u^ B (x)=2^ M (x) . (A.3) 
Eq. (|A.3j) can be rewritten as 

Xa{x)+Xb(x) = tPm(x) , (A.4) 

where 

XA (x) = (l-u/l) 1 ' 2 ^) , X b(x) = (^/2) 1/2 fcW , 

(A.5) 

which permits us to take over the analytic procedures of 
ref.[l0(. We first rotate xa(x) and xb{x) by 7r/4 in the 
function space in which they are defined, introducing 

X±(x) = ^(xa(x)± X b(x)) (A.6) 

and leaving "lengths" within that space invariant so that 

xl(x)+X 2 -(x)=i> 2 M (x) ■ (A.7) 

Finally, we introduce the polar angle [3 = (3(x) in the 
function space, 

X+(x) = tpM (x) cosP(x) , x-(x) = ipM(x)sin(3(x) , 

(A.8) 



C. The Euler equation for (3(x) 

In partition theory (PT) 0] , a Hamiltonian is assigned 
to each part for each integer number of electrons enter- 
ing into its PPLB ensemble. Since our "electrons" do not 
interact, it is sufficient to assign a one-electron Hamilto- 
nian to each part, 

2 

p 

H a = — +v a (x) , va,b(x) = -Za : b5{x ± R/2) . 

(ATI) 

The PPLB energy functional of the collection of parts is 
then 

£ = {2-v)[^ A ,H A ^ A ) + v(i> B ,H B ^ B ) 

= 2[( X a,H aXa ) + (xb,H bX b)] ■ (A.12) 

Inserting the transformation (|A.9[) and the definition 
(|A~TT|) of H a into Eq. (fA~l~2|) results in 

£ = j dx {^(x) + Vm(^) [0{xf + v A (x) + v B (x) + 

+ (v A (x)-v B (x))sm2P(x)}} (A.13) 

for the energy functional £, now a functional only of (3{x). 
In Eq. (|A.13|1 and in the following, primes indicate deriva- 
tives with respect to x. 

Varying £ with respect to /3{x) yields 

5£ = 2 J dx^ 2 M (x) {(3'(x)5(3'(x) + (v A (x) - v B (x)) 

x cos2f3(x)5(3{x)}(A.U) 
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The usual integration by parts leads to 



88 = 2{^ M {x)f3\x)5f3{x)\ x - = Z ao + 



dx 



dxV M[x) dx 



(A.15) 



+ i> 2 M {x){v A {x) - v B (x))cos2p(x)] 5(3{x)} . 
The Euler equation 

~ + ^m( x )( v a(x) - vb(x)) 

xcos2f3(x) = (A.16) 

and the boundary conditions 

ij 2 M (x)(3 , (x)5P(x)\™ oo = (A.17) 

result from imposing stationarity on £. At first glance it 
might seem that Eq. (|A.17[) is satisfied automatically since 
ip 2 M (x) J, as | oo < However, unless the boundary 
condition 



P'(x) = , |x| = oo 



(A.18) 



is imposed, (3'(x) diverges unacceptably as ip A 2 (x) as 
|z| t oo. 



The conditions (|A.20|) and (|A.21|) imply that 

ai = -Z A i; 2 M (-R/2)cos2f3 A 

= -Z B ip 2 M {R/2)cos2[3 B , \x\ < R/2(A.25) 



From Eq. (|A.20|) ct2 is (3 a when the lower limit in 
Eq. (|A.23|) is set at —R/2, yielding a second relation be- 
tween Pa and Pb, 



Pb ~ Pa = at x 

Inserting the explicit form 
gral in Eq. (|A.26j) results in 



R/2 



dx 



il> 2 M {x) 



(A.26) 



-R/2 W M \ 

for tPm{x) into the inte- 



R/2 



dx 



1 



R/2 YM 



AW 



2kD 



„kR/2 



-kR/2 



ip M (-R/2) MR/2) 

(A.27) 

Taken together, Eqs. (|A.25HA.27l) fix the values of Pa,b, 
determining p(x) for |x| > R/2, and Eq. (|A.23|) then de- 
termines P{x) for | as j < R/2. 



E. The partition potential 

Since the Xa(x) are proportional to the ip a (x), 
Eq. (|A.5[) . they satisfy the same Schrodinger Eqs. (|2. 171) . 
Summing over a and dividing by Xa(x) + Xb(x) yields 



D. Solving for (3(x) 



Because the v a {x) are 8- function potentials, Eq. (|A.16|) 
reduces to 



.. , 2 ,d(3{x) 



d_ 

dx 



dx 



= 



(A.19) 



with the additional boundary conditions 

P{-\B+)=P{-\R-) = p A \,__ip 
13'{-\R-)- l3'{-\R+) = Z A cos2p A J 2 ' 

(A.20) 



/3'(ii?-)-/3'(ii?+) = -Z B cos2p B ] x 2 



>2 " ' ^ ^2 

The general solution of Eq. (|A.19[) is 

d(3(x) ot\ 
dx iPm{x) 



(A.21) 



(A.22) 



p{x) 



F dx 'm^ + ^ ■ (A - 23) 



with the constants ai and a2 taking on different values 
in the three domains |ar| > R/2, —^R < x < ^R. The 
condition (|A.18|) implies that otx and P'{x) vanish for 
|x| > R/2 so that 

P{x) = /3 A , x < -R/2 and (3(x) = Pb , x> R/2 . 

(A.24) 



v P (x) 



Xa(x) +Xb(x) 2 
v a (x)xa{x) + v b (x)xb{x 



Xa(x) +Xb(x) 



(xa(x) +Xb(x)) + 
(A.28) 



Expressing xa{x) and xb{x) in terms of X+i x ) and 
X-( x ) via Eq. (|A.6p . and using Eq. (|A.8j) . results in 



v P {x) = flM + 



2 M x)cosP{x) dx-^ M{x)C ° Sm) + 



- ^X)"a(aO(l + s a tan/3 Q ) , (A.29) 

a 

where s a — 1 for a — A and s a = — 1 for a = B. The 8- 
function character of v a (x) and the definitions of Pa and 
Pb of Eqs. (|A.19p - (|A.20p were also taken into account in 
arriving at Eq. ([A~29|) . Since the molecular wavefunction 
ipia{x) satisfies 



1 d 2 ip M (x) 



+ (va(x) + v b (x))iPm(x) = hm^m{x) , 

(A.30) 



2 dx 2 

Eq. (|A.29[) can be transformed to 



v P (x) 



1 



i&nfiix) 
d 2 p(x) 



dipM{ x ) dp(x) 



dx 2 



^m(x) dx 
fdP(a 



dx 



\ dx 



(A.31) 



77 y^v a (x)(l - s a tan/3 Q ) 
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Using the Euler Eq. (|A.16p for (3(x), vp(x) can be further 
expressed as 

1 fd(3{x)\ 2 1 ^ . , 

x / a 

x (1 - s Q tan/3 Q (l + cos2/3 Q )) . (A.32) 

Finally, by using Eqs. (|A.23[) and (|A.25|) , we note that 
in the internuclear region \x\ < R/2, d(3(x)/dx is simply 
proportional to ip M 2 (x), yielding 

M.) = -1^ ( :5^ 2/3A W2-NIA.33) 



-^u Q (1 - s Q tan/3 Q (l + cos2/3 Q )) , 



where 8(y) (= for y < 0, 1 for y > 0) is the Heaviside 
step function. Eq. (|A.33l) correctly reduces to the parti- 
tion potential of lDi?2 [10( when Z A = Zb- 



Proof that i? OC c is proportional to InAZ for small 
AZ 



The bonding eigenvalue is 



E 



M 



-(Ea + Eb)- 



~] 2 



-(Ea-Eb) 



and 



D 
A 



Em — Ea 
Vo 



1/2 



(A.37) 



(A.38) 



Inserting Eq. (|A.37|) into Eq. l|A.38|) and rearranging gives 




(E b -Ea\ 2 


Eb-E a 


V 2« J 


2v 



(A.39) 



Now, since Eb — E A = ZAZ and «o = 
/ AZe ZR \ 2 AZe ZR 



-Z A e 



then 



B 
A 



\ 2Z 



2Z 



< 1 



(A.40) 



Thus 



We expect that i? cc, the nuclear separation at the 
crossover from covalent to ionic behavior, goes to infinity 
as AZ I 0. There, the fragment wavefunctions must 
approach the free-atom wavefunctions, and the tight- 
binding LCAO must be a good approximation to the 
molecular orbital. We can thus write 

iMs) = AiP° A (x) + Bi> B {x) (A.34) 

where the if>a( x ) are the orbitals of the isolated atoms, 
Eq. (|2.ip . Taking matrix elements of the molecular 
Hamiltonian yields equations for A, B, and the molec- 
ular energy Em- 

(A.35a) 
B = 
B = 
(A.35b) 



[E° A 
[(E° 



E 



M 



V A B A ] 



! -E m )Sab +v^ A ] A 



- Em) Sab 
\E b — Em 



ctfi 
' J A/B 



are ma- 



where the E^, are the energies of Eq. 
trix elements of the A/B potentials of Eq. (|A.ll[) , and 
Sab is the overlap {ijj A , V's)- I n evaluating all quantities 
except Ea — Em and Eb 
AZ i 0. The result is 



Em, we can take the limit 



Sab 

v AA 
V B 

V A 



Sba 

V A 
„,BA 



Sr 



ZRe~ ZR 

J _ Z 2 e -2ZR 

, -Z 2 e- ZR 



The complete solution of Eqs. (|A.35a| - (jA.35bp shows that 
S enters into Em, A, and B as 1 — S 2 and Svd which 
become exponentially small corrections as AZ J. 0, R | 
co. Similarly Vd enters only in the combination Svd and 
can be neglected as well. Thus the equations for A and 
B simplify to the classic bonding-antibonding equations 

(E A - E M )A + v Q B = (A.36a) 
v Q A+(E B -E M )B = (A.36b) 



| - f(y) < i , v i 

A = (l + /(y) 2 r 1/2 
The molecular density is 



AZ 
~2Z l 



ZB. 



n M (x) = N m iPm(x) 2 

= N M [A 2 ^ 2 A {xf + B 2 i> 2 B {x) + 2AB^ A (x)i> B {x)} 
N M Z 2 



i + f{v) 



; {gi{x) + g 2 {x)) 



(A.41) 



where 



9i(x) = e- 2Z \*\ + f\y)e- 2Z \*- R \ 
g 2 (x) = 2f(y)e~ z ^\+\*- R ^ 

It can be checked that gi[x) always exceeds g 2 (x) and 
becomes exponentially larger than g 2 (x) as x departs 
from \R - i(ln/(y))/Z. The cross term in Eq. (|A~4l|) 
can then be neglected and um(x) is thus of the form ob- 
tained from Partition Theory: tim(x) = N M A 2 ip A (x) 2 + 
NmB 2 ^P 2 3 (x), so that 2A 2 can be identified with 2 — v 
and 2B 2 can be identified with v, implying that 

f 2 (y) 



N M - v 



(A.42) 



We have chosen v — Nm/3 to define R' occ - That implies 
that f(y) = l/y/2 at R occ ' for N M = 1 (or f{y) = 1/y/E 
for Nm — 2). This leads to the observed behavior of R' occ 
in Fig.([T5]): 

R' occ = UnZ-Un{V2AZ) (A.43) 
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